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Abstract 



In this paper we study certain aspects of the complete integrability of the Genera- 
' lized Weierstrass system in the context of the Sinh-Gordon type equation. Using the 

conditional symmetry approach, we construct the Backlund transformation for the 
Generalized Weierstrass system which is determined by coupled Riccati equations. 
Next a linear spectral problem is found which is determined by nonsingular 2x2 
matrices based on an sl(2, C) representation. We derive the explicit form of the Dar- 
boux transformation for the Weierstrass system. New classes of multisoliton solutions 
of the Generalized Weierstrass system are obtained through the use of the Backlund 
\ transformation and some physical applications of these results in the area of classical 

OA ' string theory are presented. 

O 

t— i ; 

1 Introduction 

O 

The theory of surfaces has recently had a large resurgence of interest in many applications 
to very different and diverse areas of theoretical physics. This is in no small measure due 
to the fact that many equations which are of interest in various areas of mathematical 
physics also appear, or can be incorporated, in the study of surfaces in three-dimensional 
space. In particular, the Sine and Sinh-Gordon equations are of a great significance in this 
respect. Moreover, the theory of constant mean curvature surfaces has been of continued 
importance in the study of many problems which have physical applications. For example, 
minimal surfaces and constant mean curvature surfaces in particular have recently found 
applications in the areas of two-dimensional gravity [1], quantum field theory as well as 
in the area of string theory [1]. As a string propagates through space-time, it describes 
a surface in space-time which is called its world sheet. The study of classical strings in 
three and four-dimensional space is a crucial first step in producing a quantum theory of 
strings. It can be said that first quantized string theory is the study of conformal field 
theories on Riemann surfaces. 
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The stationary two-dimensional sigma model has been shown to be of use in generating 
two-dimensional surfaces immersed in three-space, and there are links between this model 
and other models, such as the non-Abelian Chern-Simons theories, which are of great 
importance in certain areas of condensed matter physics [1]. It has been shown that the 
Chern-Simons gauged Landau-Ginsburg model plays the role of effective theory for the 
Fractional Quantum Hall Effect [2] . The Chern-Simons equation of motion can describe 
time evolving two-dimensional surfaces in such a way that the deformation is not only 
locally compatible with the Gauss-Codazzi equation, but completely integrable as well [3]. 
In the static limit, the self-dual version of the model possesses soliton solutions. These 
correspond to Laughlin's quasiparticles and give a realization of anyon quasiparticles. On 
the other hand, the self-dual Chern-Simons model can be associated with the stationary 
two-dimensional continuous classical Heisenberg model, which can be related to the two- 
dimensional sigma-model. 

Another area of recent interest with regard to applications is to the area of liquid 
crystals and the theory of membranes [4]. Fluid membranes may be idealized as two- 
dimensional surfaces with each membrane being made up of a double layer of long mo- 
lecules. Various physical properties of interest such as elastic free energy per unit area 
can be calculated in terms of quantities which are directly related to the geometry of the 
surface. In fact the curvature elastic free energy per unit area of the membrane can also be 
formulated rigorously in terms of two-dimensional differential invariants of the surface. It is 
of considerable interest with regard to these types of physical applications to obtain shape 
equations for the membrane. These interrelate the basic parameters and functions which 
determine the form of a given membrane surface, as in a liquid crystal. In an equilibrium 
state, the energy of any physical system must be minimized. One usually writes down 
a shape energy function F in terms of the basic parameters and then minimizes it, and 
the result is a shape equation. An example of such a shape function is given by 



where k c is the bending rigidity of the membrane, H the mean curvature and the spon- 
taneous curvature Co takes account of the asymmetry effect of the membrane or the sur- 
rounding environment. The pressure difference between the outside and the inside of the 
membrane is called Ap, A the tensile stress acting on the membrane. Mathematically, Ap 
and A may be considered as Lagrange multipliers. The shape equation is obtained from 
the first variation of this F. Specific Delaunay's surfaces of constant mean curvature can 
be written as sin^(p) = ap + d/p. This equation is then substituted into a given shape 
equation and results in constraint equations between the parameters a and b, which give 
a characterization of the surface [4, p. 114]. 

In this paper, we study a connection between the Generalized Weierstrass (GW) sys- 
tem inducing constant mean curvature surfaces immersed in M 3 and a Sinh-Gordon type 
equation. The objective of this paper is, using this link between these two systems, to 
derive a representation of a linear spectral problem for which the matrices are nonsingular, 
and find the corresponding Darboux and Backhand transformations for the GW system. 
These transformations are derived here for the first time. Based on these transformations, 
we construct solutions of the GW system and investigate minimal surfaces immersed in 
three-dimensional Euclidean space. 
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This paper is organized as follows. In Section 2, a short presentation of the conditional 
symmetry approach for partial differential equations which admit the Painleve property 
is given. In Section 3, using the conditional symmetries we derive the linear spectral 
problem, and we find the Darboux and Backlund transformations for the GW system. 
Section 4 contains new examples of multi-soliton solutions of the GW system and some 
physical interpretations of these results are given in the area of classical string theory. 

2 Conditional symmetries 

In this section, we give a brief overview of the conditional symmetry approach for PDEs 
as developed in [5-8]. In this context, we concentrate on examining certain aspects of 
integrability of fc-th order nonlinear PDEs. The technique which is outlined below is 
applied only to such classes of PDE which pass the Painleve test and can be presented 
in the form of a polynomial in the unknown variable u and its derivatives, possibly after 
a transformation in the space of independent and dependent variables X x U. The basic 
terminology and notation used here in the application of Lie groups to differential equations 
are in conformity with [9] . We are particularly interested in combining singularity structure 
analysis and Lie point symmetries in order to recover the Auto-Backlund transformation 
(Auto-BT) and Darboux transformation (DT) for PDEs if such exist. In the literature, 
several attempts at treating this subject can be found, for example [5-8], and references 
therein. More recently, this subject has been studied for first order systems of PDEs, by one 
of the authors, leading to the development of a new version of the conditional symmetry 
method [10, 11]. This approach, like other nonclassical methods (see for a review of 
the subject [12]) makes possible the explicit determination of certain classes of solutions, 
invariant under a group of transformations which maps a subset of solutions of the initial 
equation into other solutions of a different equation, that is, a subsystem composed of 
the initial PDE and differential constraints (DCs), which are mutually consistent. In this 
presentation, we postulate in accordance with the method worked out in [11], that these 
multiple DCs take a specific form for which all derivatives of the unknown function u are 
expressible in terms of some functions of the independent and dependent variables only. 
Hence we consider the overdetermined system composed of a nondegenerate fc-th order 
scalar PDE and a first order system of DCs 



The brackets [i,j] denote the alternation with respect to the indices i and j, that is, 

4>[i,j] = 2(4>i,j - <f>j,i), 
<f>[i<f>j],u = 2(<Pi4>j,u - 4>j<Pi,u)- 




(2.1) 



(2.2) 
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Note the equation (2.2) means that the characteristics Qi of a set of p-linearly independent 
vector fields (defined on X x U) 

Zi = d x i + (f)i(x,u)d u , i = l,...,p, (2.4) 

are equal to zero. 

An Abelian Lie algebra L spanned by the vector fields Z\ , . . . , Z p is called a conditional 
symmetry algebra of the k-th order PDE (2.1), if the vector fields Z±, . . . ,Z p are tangent 
to the subvariety 

s = 5a n Sq, 

where we associate the initial system A : J k — > R and a first order system of DCs 
Qi : J 1 — > R p with the subvarieties of the solution spaces 

Sa = {(x,u^) G J k : a(x,u^) = o}, 
S q = { (x, «W) £ J 1 : Qi (x, «W) = 0, i = 1, . . . ,p} , 

respectively. This definition means that the fc-th prolongation of the vector fields Zi 
belongs to the tangent space to S at [x,u^), that is, 

prWZils € T (X)tt(fc)) 5, i = l,...,p. (2.5) 

A solution u = /(x) of the k-th order PDE (2.1) is called a conditionally invariant 
solution, if its graph {(x,/(x))} is invariant under an Abelian distribution of the vector 
fields Z±, . . . , Z p satisfying conditions (2.5). 

It has been shown [11] that a nondegenerate k-th order PDE (2.1) admits a p-dimensio- 
nal conditional symmetry algebra L if and only if there exists a set of p linearly independent 
vector fields (2.4) for which the C k ~ l functions (pi satisfy the conditions (2.3). The graph 
of a solution of the overdetermined system composed of (2.1) and (2.2) is invariant under 
the vector fields Zi, 1 < i < p. Hence according to the above definition, this means that 
there exists a conditionally invariant solution of PDE (2.1). 

It has been proved [13] that any PDE (2.1) of the fc-th order admits infinitely many 
compatible first order DCs. However, this statement shows only existence of such con- 
straints, but does not provide a constructive method for finding the explicit form of these 
DCs. Thus, the construction of conditional symmetries is reduced to the selection of such 
subsystems composed of initial PDE (2.1) and DCs (2.2) for which conditions (2.3) hold. 
In general, system (2.3) is a nonlinear one and usually very difficult to solve, except in some 
particular cases. Nevertheless, there exist many physically interesting systems of PDEs 
for which particular solutions of (2.3) lead to Backlund transformations described by first 
order differential constraints or to solutions depending on some arbitrary constants [14]. 
These particular solutions of (2.3) are obtained by expanding each function fa into a poly- 
nomial in the dependent variables u. This polynomial is reduced often to a trinomial. 
This means that equations (2.2) become a Riccati system of equations which possess the 
Painleve property. 

In Section 3, we show on a specific example of the generalized Weierstrass system 
inducing constant mean curvature surfaces in R 3 , that the conditions (2.3) for the existence 
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of conditional symmetries can be used to construct a certain class of Backlund and Darboux 
transformations. This construction consists of the following steps. 

First, we assume that PDE (2.1) passes the Painleve test. This means that (2.1) 
satisfies the necessary conditions for the absence of movable critical singularities in its 
general solution on arbitrary noncharacteristic surfaces [15, 16]. This test provides us 
with a tool for assessing the integrability of PDEs. At least we can reject as nonintegrable 
these PDEs which do not pass the Painleve test. 

We restrict our considerations to the particular case when the singularity structure 
of a solution of PDE (2.1) consists of only poles. According to Painleve analysis for 
PDEs [16], it has been shown that if the choice of the expansion variable x m terms of the 
singular manifold variable <p(x) = fo, (where x vanishes as ip — <po), then the beginning 
of the Laurent series for the solution u takes place at the order — n in that is, 

u ~ X ~ n , n € Z < 0, (2.6) 

where the power (— n) denotes the multiplicity of a pole. 

Secondly, we postulate that the difference of two distinct solutions u and u of PDE (2.1) 
can be represented in a polynomial form in terms of an auxiliary variable y up to the 
degree (— n). The variable y is a mapping of space X into some m-dimensional space B. 
This demand implies that we perform an embedding transformation of the variable (u — u) 
into the space B with coordinates y = (y 1 , . . . , y m ) . Thus we assume that it can be realized 
through a specific Darboux transformation of the form 

u-u = J2cjy J , l<#J<(-n), (2.7) 
J 

where J = (ji, ■ ■ ■ ,j m ) and ji is a positive integer such that 

#J = ]! + ■■■+ j m = ~n. (2.8) 

The coefficients c,j in the expansion (2.7) are assumed to be constants. Eliminating the 
function u in the initial equation (2.1) through transformation (2.7), we obtain a k-th 
order PDE for the unknown function y. Denote this PDE by 

A 1 (x,uW,yW)=0, (2.9) 

where vP^ is a given function of x. Equation (2.9) is the starting point for our further 
analysis. In this section, we focus on a close connection between the conditional symmetries 
and Backlund transformations associated with nonlinear systems of PDEs (2.9). This 
connection is mainly due to the fact that the set of DCs (2.2) in the variable y admits 
a superposition formula (SF), as is also the case for BTs [17]. Based on the application of 
Lie's theorem on fundamental sets of solutions [18] to the case of first order PDEs [17], we 
show this link on a specific example in the next section. We assume that these DCs take 
the particular form of coupled matrix Riccati equations based on some given representation 
of the Lie algebra. This means that the DCs (2.2) take a specific form for which the first 
derivatives of y are decomposable in terms of x and y as follows 

^l( x ) = j2A\(x)bf(y{x)). (2.10) 



234 



P Bracken and A M Grundland 



The set of functions bf is identified with linearly independent vector fields on a space B 

bi = bf(y)d ya , l = l,...r, (2.11) 

which generate a finite-dimensional Lie algebra Q 

[k,h] = C? k h, l<a,k,l<r, (2.12) 

where Cf k are constants of the assumed Lie algebra structure. Note that the problem of 
construction and classification of all finite-dimensional Lie algebras which can be realized 
in terms of operators (2.11) remains an open one [18-20]. Nevertheless, this subject has 
been recently extensively investigated and there exist numerous lists of finite-dimensional 
Lie algebras which can serve as a source for our selection (see for example [19, 20] and 
references therein). 

We select one of the finite-dimensional Lie algebras Q and its representation in terms of 
vector fields (2.11) with polynomial coefficients in the variable y. We start with the lowest- 
dimensional algebra and, if necessary, proceed to consider the higher dimensional ones. 
For a chosen Lie algebra Q the right hand side of expression (2.10) becomes a polynomial in 
the dependent variable y. Substituting (2.10) repeatedly (k — 1) times into PDE (2.9) and 
next requiring that the coefficients of successive powers of y in the equation so obtained 
vanish, we get a system of (k — 1) order PDEs for the functions A\. Denote this system 
by 

Q a (xM k) {x),Af~ l) ) =0, a = l,...,r. (2.13) 

For the assumed Lie algebra Q the compatibility conditions for DCs (2.10) impose zero 
gauge curvature conditions on the functions A\ 

A [i,J] + \ C lk A \ A ) = 0, i^j = l,...,p. (2.14) 

As a result, we arrive at an overdetermined system of equations for the functions A\. 
We denote this system, consisting of equations (2.13) and (2.14), by H. To establish 
the existence of solutions of this overdetermined system, an analysis of the compatibility 
conditions of the system is required. 

It has been shown [21] on several examples, such as the AKNS, Boussineq, Kadomtsev- 
Petviashvili, Sawada-Kotera, and Tzitzeica equations, that there exist many physically 
interesting systems of PDEs for which the general solution of the system H lead to Auto- 
Backlund transformations. This phenomena takes place when the set of solutions of the 
system H is parametrized by a function u satisfying the original PDE (2.1) and by at 
least one constant parameter A. In this case, the functions A\ can be expressed uniquely 
in terms of the old solution u of (2.1) and some constant A. Then DCs (2.10) become 

^{x) = j^A\{u,\)bt{y{x)). (2.15) 
i=i 

From each solution u of PDE (2.1) we integrate DCs (2.15) for the functions y a and 
find a solution u of original equation (2.1) via Darboux transformation (2.7). Thus, sys- 
tem (2.15) together with (2.7) determines a specific Auto-BT between sets of solutions of 
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the initial PDE (2.1). By eliminating the auxilliary variables y a from the system (2.15) 
and (2.7) we can obtain an explicit form of this Auto-BT. 

Note that the solution of (2.9) which is obtained from the proposed procedure con- 
stitutes a conditionally invariant solution since it is invariant under the p-dimensional 
conditional symmetry algebra L spanned by the vector fields of the form 

r 

Zi = d xi +J2 4(^,^1, i = l,...,p, (2.16) 
i=i 

where b\ = b^(y)d y a generate a finite-dimensional Lie algebra (2.12). 



3 The generalized Weierstrass system 

We now proceed to apply the conditional symmetry approach to the case of the generalized 
Weierstrass system (GW) inducing constant mean curvature surfaces embedded in R 3 as 
derived by B Konopelchenko in [22]. This system is described by a set of Dirac type 
equations for two complex fields ipi and ■02 given by 

dtp! = pip 2 , dtp 2 = -ptpi, 

dip! = ptp 2 , dtjj 2 = -P0i, (3.1) 

where p = |0i| 2 + \ip2\ 2 , 9 = d/dz and 3 = d/dz. 

Equations (3.1) possess several conserved quantities 

d (0i dtjj2 - ip2dipi) = 0, d (tpidip2 - fadipi) = 0, 

a(0!) 2 + a(0 2 ) 2 = o, 8 (0i) 2 + d (0 2 ) 2 = o. (3.2) 

It has been shown [22] as a consequence of these conservation laws that there exist three 
real valued functions X l (z, z), i = 1, 2, 3 such that 

X 1 + iX 2 = 2i j (0 2 dz' - i?2 dz') , X 1 - iX 2 = 2i J dz' - ^\ dz) , 

X 3 = -2^ (0!0 2 dz' + Vi02 dz') , (3.3) 

Note that by virtue of the conservation laws (3.2), the right hand sides in expression (3.3) 
do not depend on the choice of contour T in the complex plane C, but only on its endpoints. 
This is due to the fact that the integrals (3.3) have the form 

F(z, z) dz + F(z, z) dz, 

which satisfy the condition 

dF = dF. 

Consequently, the differentials of equations (3.3) are exact ones. The functions X % (z,z) 
can be identified as the coordinates of the position vector J of a surface embedded in R 3 . 
Hence from equations (3.1) and (3.3), we can determine the induced metric of this surface 



ds 2 = 4p 2 dz dz, 



(3.4) 
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in isothermic coordinates. The Gaussian curvature and constant mean curvature can be 
evaluated from K = —p~ 2 dd(lnp), H = 1, respectively, where p is given in (3.1). 

The system (3.1) is known to be completely integrable, and a Lax pair for it has been 
found recently by authors [23, 24]. However, the matrix appearing in the Lax pair is sin- 
gular and nilpotent, which prevents the construction of solutions of (3.1). The objective of 
this paper is to demonstrate that it is possible to find such a representation of the Lax pair 
for which the matrices are nonsingular. This representation is suitable for constructing 
solutions of GW system (3.1). Next, we construct the Backhand transformation and the 
Darboux transformation by making use of a close connection between these transforma- 
tions and conditional symmetries. Finally, based on these transformations, we construct 
several new classes of multisoliton solutions for (3.1) and investigate the corresponding 
surfaces in M 3 . 

In our investigation of the integrability of the GW system (3.1) we change the dependent 
variables ipi and ip2 to new dependent variables 

p= |^i| 2 + |^ 2 | 2 , 

and the current 

J = ipidfo - fodipi, (3.5) 

in order to simplify its structure. In terms of these new variables, we show that GW 
system (3.1) can be decoupled into a direct sum of the elliptic Sh-Gordon type equation 
and the conservation of current dJ = 0. 

In fact, differentiating the function p with respect to z and z, we obtain, 

dp = ipi (dtpi) + tp 2 (di/j 2 ) , Bp = tpi (&if)i) + ip 2 (9^2) , 

ddp = dipidipi + 8^29^2 — p 3 - (3.6) 

Making use of conservation laws (3.2) and equations (3.6), GW system (3.1) takes the 
decoupled form in the variables p and J, 



ddlnp= i-i- -p 2 , 8J = 0, 8J = 0. (3.7) 



If we introduce the new dependent variable 

p = e^ 2 , 

into equation (3.7), we then obtain an elliptic sinh-Gordon type equation of the form [25] 

ddip = -4sinhv2-2(l- | J| 2 ) e~^, dJ = 0. (3.8) 

In particular, if the modulus of the current J is different from zero, \J\ 7^ 0, then we can 
introduce new independent, dependent variables rj, fj and u 

dij = J 1 / 2 dz, dfj = J 1 / 2 dz, uj = — -, 

respectively, such that GW system (2.1) takes the decoupled form 

{]nu)rfi = 2 - wj , d,J = 0, dJ = 0. 



On Complete Integrability of the Generalized Weierstrass System 



237 



Hence GW system (2.1) becomes a direct sum of the elliptic Sh-Gordon equation and the 
conservation of current J. 

Equation (3.7) has the Painleve property. Consequently, we obtain that the general 
solution p of (3.7) admits double poles with two residues of opposite sign 

P ±2 = = ±x~\ (3-9) 

where \ is the expansion variable of the Laurent series. According to the proposed pro- 
cedure, based on the approach described by Conte and Musette [26] we assume a specific 
form of the Darboux transformation for the case when equation (3.8) admits opposite 
residues 

ip-v = 2\n— . (3.10) 

The function v satisfies the initial equation (3.8) and quantities <j)\ and (f>2 are two entire 
functions. Note that a similar situation for double poles with opposite residues arises in 
the study of singularity structure for sine-Gordon and MKDV equations [26] . Introducing 
a new variable y = (fii/4>2 and changing the variables in (3.10) according to 

p = e^ 2 , q = e vl \ (3.11) 

we find that the Darboux transformation for equation (3.7) can be realized by the following 
expression 

p = qy, (3.12) 

or equivalently, 

ip = v — 2 In y. 

A first step on the way to constructing a BT for (3.7) is to look for a conditional symmetry 
algebra L spanned by two vector fields Z\, Z 2 which have the characteristic equations of 
the form (2.10). We have to assume a specific Lie algebra structure for the generators {6/}. 
We start the analysis with the lowest-dimensional case, namely that of the sl(2, C) algebra 
which admits the one-dimensional representation (B y , yB y , y 2 B y ) in terms of a coordinate y. 
This algebra comes up in the study of several completely integrable models eg. [19]. In 
our case, DCs (2.10) in one complex variable y take the form of coupled scalar Riccati 
equations with nonconstant coefficients 

dy = A\(z, z) + A\(z, z)y + A 2 (z, z)y 2 , 

By = A° 2 (z, z) + Al(z, z)y + A 2 (z, z)y 2 . (3.13) 

The zero curvature conditions for (3.13) are given by 

3A\ - 8A° 2 + A\A° 2 - A\A\ = 0, 
BA\-dAl + 2(A 2 A 2 -A 2 A° 1 )=0, 
BA\ - dA\ - A\A\ + A\A\ = 0. (3.14) 
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The substitution of the new variables (3.11) and the Ansatz (3.12) into equations (3.7) 
gives 

— (ddq y + BqBy + BqBy + qddy) \- (dqy + q dy) (Bq) 

qy q 2 y 

1 - Ul 2 

-— (Bqy + qBy)By-^ + q 2 y 2 = 0. (3.15) 



Using equations (3.13) we can eliminate the derivatives of the complex variable y in the 
expression (3.15). Next, we require that the coefficients of the successive powers of y in 
the equation so obtained vanish, to give the system 

(1) q 2 A° 1 A°2 + \J\ 2 = 0, q 2 A 1 A°2 + \J\ 2 =0, 

(2) BA° 2 - A\A° 2 = 0, BAl - A^Al = 0, 

(3) dA 1 2 + AlA° 1 -A 2 1 Al + dd\nq = 0, BA\ + A\A\ - A\A\ + B8\nq = 0, (3.16) 

(4) dA\ + A\A\ = 0, BA\ + A\A\ = 0, 

(5) A 2 A 2 + q 2 = 0, A 2 A 2 + q 2 = 0. 

We obtain an overdetermined system (denoted by H in the previous section) composed 
of (3.14) and (3.16) for the unknown functions A\. In our case, this system is consistent 
since the compatibility conditions are identically satisfied. This system has a nontrivial, 
unique solution for the A\. We briefly outline how this solution can be obtained. The 
second and fourth equations in the first column of (3.16) can be written in the form 

B\nA^ = A\, -<91n^ = ^. 

Equating these equations, we can integrate to obtain 

Al = (3.17) 

where g is a complex function of z. Similarly, from the second and fourth equations in the 
second column of (3.16), we obtain 



which can be solved to give, 



9 In A? = -din A?, 



A? = ^, (3.18) 



where h is a complex function of z. Substituting these results into the first equation 
in (3.16), we find that 

h(z)g(z) - \J\ 2 =0. 

Thus, without loss of generality, one may take h(z) = J(z) and g(z) = J(z). 
From the first equation in (3.16), we can write 
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thus A\ is determined in terms of A®. Equating equations (3.18) and (3.19), we elimi- 
nate A® and we get, 

A\ = -t A l (3.20) 

From this, we can substitute A\ into the fifth pair of equations in (3.16), and we 
obtain A\ in terms of A\, as 

A\ = j,- (3.21) 

Using the first column of (3.16), we substitute the relation A\ = din A® from equation 
(3.16-2) as well as (3.19) through (3.21) into equation (3.16-3). In this way, we can 
eliminate the coefficients A\, A\, A\ and A\ from differential equation (3.16-3). We 
obtain a partial differential equation for the function A\ in the form 

ddln A° - * {Alf + -g^ + ^ - <? = 0. (3.22) 
■J q 2 (A%) Q 

Moreover, if we introduce a new variable defined by 

Q = $, (3.23) 
into (3.22), it is transformed into the following form 

I /I 2 

3dlnQ- l -^ + Q 2 = 0. (3.24) 

This coincides with equation (3.7). Substituting (3.19)-(3.21) into the pair of Riccati 
equations (3.13), we obtain 

dy = -^ + d\nA°y- q ^3y\ B.J = 0, 

dy = Al + B\nl^^jy+^y\ 8 J = 0. (3.25) 

The compatibility condition for (3.25) is satisfied identically whenever (3.7) holds. 

Note that Q is some solution to (3.7) which is related to q by (3.23). We could obtain 
a particular form for ^4° by considering the case in which Q = q. Then (3.23) implies that 



A° 2 = iVJ. 

Equation (3.7) is invariant under the transformation 

J -► AJ, (3.26) 
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provided that |A| 2 = 1. A Backlund parameter A can be introduced into equations (3.25) 
by carrying out transformation (3.26). In this case, the pair of Riccati equations (3.25) 
takes the form 



By 




By = fX 1 / 2 JV2 + Bin 



'jl/2 



8J = 0, 



y-zAVV^V, |A| 2 = 1, 



(3.27) 



where q satisfies (3.7). Hence the above DCs become an Auto-BT for GW system (3.7), 
while the Darboux transformation is defined by (3.12). Furthermore, by linearizing the 
Riccati system (3.27), we obtain the associated linear spectral problem for (3.7) with 
spectral parameter p, of the form 











d 



( 

V - v./. 

/ i- j 1 / 2 



ip 



1/2 



J 1 ' 2 J \ 



iq< 



(Sf 



V 



m 
l 



1/2J1/2 



71/2 



) 



(3.28) 



where y = 4>i/4>2, |/u| 2 = 1 and Bj = 0. The Lax pair (3.28) is based on a nondegener- 
ate sl(2,C) representation. Note that for any holomorphic function J, the compatibility 
condition for (3.28) reproduces the system (3.7) in the variable q. 



4 Multi-soliton solutions 

A number of new types of soliton solutions to the GW system (3.1) will now be presented 
here based on the BT (3.27). Let us define a new dependent variable 

P=f (4.1) 

It has been shown in [23], that if functions ip\ and ip2 are solutions of GW system (3.1), 
then the function p defined by (4.1) is a solution of the two-dimensional Euclidean sigma 
model equations 



ddp--^-^dpBp = 0, d Bp--^-^dpBp = 0. (4.2) 
1 + \p\ z 1 + \p\ z 



Conversely, if p is a solution of the sigma model (4.2), then the solutions V'l and ^2 of 
GW system (3.1) have the form 

(Bp) 1 / 2 (dp) 1 / 2 
^ = e ^iZT~|2 ' ^ = ef^, 6 = ±1. (4.3) 
1 + \p\ z 1 + \p\ z 
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Thus solutions to GW system (3.1) can be obtained directly by applying the transforma- 
tion (4.3) when a solution of the sigma model (4.2) is known. Once we have found particu- 
lar solutions i/ji, we can calculate the corresponding value for the function q = IV'i | 2 + IV^ | 2 - 
Next, we employ succesively the Auto-BT (3.27) in order to find new solutions p and sub- 
sequently from (3.1), we construct the multi-soliton solution ipi for GW system (3.1). 

1. First, we look for simple nonsplitting rational soliton solutions of (4.2) admitting 
one simple pole given by [27] 

Pj = - — r 1 , j = 1, • • • , N, aj £ C. 

Z-ttj 

Using the Auto-BT (3.27), we get the following algebraic multi-soliton solution of (3.1) 

■■■■■ :>:'... n 



^ (z — a,-) z — au 
1=1 ^ ■>' k=i K 



N , N 



j=i y j ' k=i 



±i. 



For N = 1, the surface is determined by the equation 

The corresponding surface of revolution with constant mean curvature H = 1 for a\ = 2 
is plotted in Fig. 1. Such a surface in R 3 with a similar shape has recently been obtained 
in a cosmological application to white hole fissioning [28]. 

2. A large class of hyperbolic nonsplitting solutions, (ddp ^ 0), of the sigma model 
equations (4.2) can be constructed when the function p satisfies the algebraic constraint 
\p\ 2 = 1. Consider a class of nonsplitting hyperbolic solutions of (4.2) 

N 



p = exp(cosh(z — cii) — cosh(z — a^)). 



i=i 



In all of the solutions of (4.2) which are presented below, the Oj will be arbitrary complex 
constants. The derivatives of p with respect to d and 8 are given by 



N N 



dp = ^ sinh(z — aj)p, dp = — ^ sinh(z — a,j)p. (4.5) 



i=\ i=l 



respectively. Substituting (4.5) into (4.3), we obtain the following solutions to GW sys- 
tem (3.1), 

(N \ 1 / 2 / N \ 1 / 2 

p^sinh(z - oi) J , i/j 2 = g I P^sinh(z - at) J 

i N 

^sinh(z-ai) . (4.6) 



P= 2 



i=i 
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Figure 1. Constant mean curvature surface corresponding to equation (4.4) for a\ = 2. 

Note that these solutions do not admit any singularities. Now, solution (4.6) for N = 1 is 
substituted into integrals (3.3), we obtain 



Eliminating the z-dependent factors on the right of (4.7), the following relationship be- 
tween the Xi variables is obtained 



This represents a sphere of unit radius. Note that similar results hold when sinh is used 
in place of cosh in expression (4.4). 

3. Consider a class of hyperbolic nonsplitting solution of (4.2) which are obtained from 
the tanh function that satisfies the algebraic condition \p\ 2 = 1. This type of solution 
represents kink-type solution, and is generated by 



X\ + 1X2 = i sinh(— cosh(z — a) + cosh(z — a)), 
X\ — iX<i = — i sinh(— cosh(z — a) + cosh(z — a)), 
X3 = — cosh(cosh(z — a) — cosh(z — a)). 



(4.7) 



x\ + x\ + x\ = 1. 



(4.8) 



A? 




(4.9) 



i=i 



The derivatives of p in (4.9) are 



N 



N 





(4.10) 



i=i 



i=i 
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f N \V2 

p^sech 2 (z - at) I 

V 1=1 / 



Substituting (4.10) into (4.3), we obtain the following multi-soliton solutions tpi to GW 
system (3.1), 

(N \ 1 / 2 / N \ 1 / 2 

pJ2 sech 2 (z - ai)\ , ^2 = 

1 N 

p= 2 J2 secl]2 ( z - a i) ■ ( 4U ) 
i=\ 

Note that the functions V« admit only simple poles. For N = 1 the surface represents 
a sphere of radius one. 

4. All functions p which generate solutions of (4.2) in Examples 1 and 2 can be used to 
generate larger classes of solution by taking the functions which appear in the sums for p 
in expressions (4.4) and (4.9) and combining them by taking products in different ways. 
For example, there exists a hyperbolic, nonsplitting solution of (4.2) of the form 



/ N \ 

p = exp y~^(cosh(2: — a,i) — cosh(z — ai) + sinh(z — a«) — sinh(z — m) . 
The corresponding solution of GW system (3.1) has the following form 

(N \ 1 / 2 

p^2(s'mh(z - ai) + cosh(z - ai)) \ 

/ JV \ V2 

^2 = - f /9^(sinh(z - ai) + cosh(z - aj))J 



(4.12) 



N 



y~^(sinh(z — Oj) + cosh(2; — a«)) 



i=i 



(4.13) 



This solution does not admit any singularities. For N = 1, the surface represents a sphere 
of radius one exactly of the form (4.8). 

5. Let the complex function gi(z,z) be a set of N harmonic functions and fi(z) a set 
of N arbitrary complex valued functions of one complex variable z, and the complex 
conjugates of these. Then we can write a general solution of the sigma model as follows 



exp 



N 
3=1 



N 

n 



Making use of transformation (4.3), the corresponding general solution of GW system (3.1) 
has the form 

1/2 



N 



N 



i>2 = 



N 



i=i 3 

N 



dfj(z) 

\ m 

1/2 



E , \- dfj(z) 



±1. 



(4.14) 
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6. As a final example, let us consider a trigonometric solution generated by p of the 
form 

p = \sm(nz), n € Z, A G C. 
Using (4.3), the corresponding periodic solution of GW system (3.1) has the form 

. . (nX cos(nz)) 1 / 2 , (raAcos(nz)) 1 / 2 

ipi = e\sm(nz)- „ . — , V2 = e- 



1 + |A| 2 sin(nz) sin(nz) ' " 1 + |A| 2 sin(nz) sin(nz) 

The corresponding constant mean curvature surface for A = 1 is given by 

Y-2 , v2 _ A3 (A3 - 2) 2 
Xl+X2 " 4- A3 ' 

which has the shape of a 'Mexican hat'. It is obtained by rotation of the curve 

A 2 (4-A 3 )-X 3 (A3-2) 2 = 0, 

for A3 € [0,4), around the Ai or A2 axis. 

It is useful at this point to mention some physical applications of some of these surfaces. 
Cylinders and spheres have applications to certain types of cosmological models, and 
should also be useful in describing event horizons in general relativity [1]. In string theory, 
the motion of a particle is described by a surface which propagates through space-time. 
Minimal surfaces can constitute a way of describing particle states. 

Certain types of soliton solutions which are localized in time as well as in space are 
referred to as instantons. Such solutions exist in gauge theories, since the gauge-field 
equations are relativistic and give rise to topological nontriviality in time as well as in 
space [29]. The application of the generalized Weierstrass system to strings in three- 
dimensional Euclidean space will be outlined. This system allows the construction of any 
surface in R 3 , where p is a real- valued function of z and z. With u = \il>i\ 2 + IV^I 2 ; the 
Gaussian curvature is K = —u~ 2 dd(\ogu) and the mean curvature is H = p/u, so that 
when the mean curvature is constant, p = u. In terms of the variables p, i\>i the required 
action for the string has the following form 

S = 4p J (\ipi\ 2 + \i>2\ 2 ) 2 dxdy + ^- j p 2 dxdy. 

Classical configurations of strings can be described by common solutions of this Nambu- 
Goto-Polyakov action S and the generalized Weierstrass equations, which provide surfaces 
in R 3 . In generic coordinates, the corresponding Euler-Lagrange equation has the form 

AH + 2H (H 2 - K) - 2a fi H = 0, (4.15) 

where A is the Laplace-Beltrami operator, and under the conformal metric 

AH = u~ 2 ddH. 

In terms of the variable tp = H' 1 , p = u/ip, the Euler-Lagrange equation takes the form 



ddf + [2p 2 + ddlnp 2 } ip - 2a /«opV = 0. 
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When the mean curvature is constant, ip = tpo, and the Euler-Lagrange equation reduces 
to a second order linear equation 

dBlnp 2 + 2(l-a oW ^)p 2 = 0, 

which can be transformed into the Liouville equation 

dde + Pe e = 0, (4.16) 

where 6 = lnp 2 . For (3^0, the solution of (4.15) has the form 

fl = lndG' + ln<9G-21n(|G| 2 + /?/2) , 

where G{z) is an arbitrary analytic function. Solutions of (4.14) which correspond to mean 
curvature different from zero can be invesigated as well. This task will be undertaken in 
a future work. 

5 Final remarks 

In this paper, we demonstrate that the task of finding large classes of solutions of PDEs is 
related to the group properties of an overdetermined system composed of an initial system 
of PDEs (2.1) subjected to DCs (2.2). We show that this has a group theoretical interpreta- 
tion in terms of "conditional symmetries" . The main difficulty in this approach is related 
to finding reasonable ansatzes that yield compatible solutions. The approach adopted 
here is based on multiple constraints satisfying several specific conditions, such as (2.10), 
(2.12), (2.13) and (2.14), which enables us to overcome these difficulties. The proposed 
approach simplifies the task of solving the nonlinear determining equations (2.13) for the 
initial system (2.9) by providing us with an almost entirely algorithmic procedure. The 
most important advantage of this method is that it gives us effective tools for constructing 
certain classes of BT described by first order differential equations in a systematic way. 
Its effectiveness has been demonstrated by the results, both reconstructed [10, 21, 24, 25] 
and new in Section 3 and in Section 4 for the GW system (3.1). 

In conclusion, we have learned from this example that it is useful to subject system (2.9) 
with another one, which involves certain auxiliary variables dictated by the possible rep- 
resentations of the "hidden" symmetry algebra (2.12). Selection of appropriate algebras 
and their representations is still an open problem, requiring further investigation. Finally, 
the presented example of GW system (3.1) suggests that there may exist a connection 
between the conditional symmetry method and the prolongation structure approach of 
Wahlquist and Estabrook [30]. Indeed, the auxiliary function y is nothing else than the 
pseudopotential introduced in their approach. 
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